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自 适应 滤波 结合 改进 T-Snake 模型 的 甲状 腺 超声 图 自动 分 割 方法 


周 春 瑜 ， 程 显 煞 
(南通 大 学 电气 工程 学 院 , 江苏 南通 226006) 


摘 要 : 提出 一 种 基于 T-Snake 模型 的 甲状 腺 超声 波 图 像 分 割 的 新 方法 。 首 先 ， 结 合 基于 窗口 的 各 向 异性 扩散 滤波 
方法 与 自 适应 加 权 中 值 滤波 算法 有 效 地 消除 甲状 腺 超声 波 图 像 斑 点 噪声 ;其 次 ， 以 传统 TSnake 模型 为 基础 ， 增 加 自 
适应 的 区 域 能 量 和 膨胀 力 对 非 连 续 边 界 与 弱 边界 进行 有 效 提 取 ， 实 现 甲状 腺 超声 波 图 像 的 自动 分 割 ; 最 后 ， 设 定 模 
型 参数 ， 使 用 临床 数据 进行 实验 。 结 果 证 明 ， 应 用 该 方法 得 到 的 自动 分 割 结果 的 平均 相对 差异 度 小 于 5%， 平 均 相 
对 重合 度 大 于 91%， 验 证 了 方法 的 可 行 性 。 
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Automatic segmentation of thyroid ultrasound images based on 
adaptive filtering combined with T-Snake model 


Zhou Chunyu, Cheng Xianyi 
(School of Electrical Engineering, Nantong University, Nantong Jiangsu 226019, China) 


Abstract: This paper proposed a new segmentation method of thyroid ultrasonic image based on T-Snake model. Firstly, a 
window-based anisotropic diffusion filter, combined with an adaptive weighted median filter, effectively eliminated speckle 
noise in thyroid ultrasonic images. Secondly, the method, based on the traditional T-Snake model and the adaptive region 
energy and expansion force, extracted the discontinuous boundary and weak boundary effectively, and realized the 
automatic segmentation of thyroid ultrasonic image. Finally, it set the model parameters and experimented with the clinical 
data. The results showed that the average relative difference is less than 5% and the average relative overlap is more than 
91%. It verifies the feasibility of the proposed method. 

Key words: thyroid ultrasonic image; image segmentation; t-snake; adaptive; weighted median filter; anisotropic diffusion 


filter 
0 引言 散 滤波 方法 有 效 去 除 斑点 噪声 ;其 次 ,在 传统 T-Snake 模型 中 
加 引入 自 适 应 的 区 域 能 量 和 膨胀 力 来 提升 方法 的 鲁 棒 性 ， 实 现 
医学 图 像 的 分 割 问题 一 直 是 难题 ， 也 是 研究 热点 。 目 前 非 连续 边界 与 弱 边界 的 有 效 提取 。 


的 分 割 方法 的 精度 难以 满足 应 用 要 求 。 由 Kass 等 人 M1 提出 的 


1 ”斑点 噪声 自 适 应 的 滤波 


全 ”Snake 模型， 是 近年 来 图 像 处 理 中 用 于 目标 轮廓 提取 的 一 种 

一 主流 模型 ,但 是 ， 传 统 Snake 方法 存在 许多 不 足 ， 其 分 割 结果 原始 甲状 腺 超声 波 图 像 具 有 显著 的 斑点 噪声 ， 对 于 图 像 
对 初始 位 置 依赖 较 大 ， 对 拓扑 结构 复杂 多 变 的 轮廓 线 提取 效 ”分割 影响 较 大 ， 因 此 必须 对 原始 图 像 进 行 预 处 理 。 各 向 异性 
果 较 差 。 于 是 有 许多 研究 人 员 对 Snake 模型 进行 了 改进 ?> 7， ”扩散 算法 (SRAD) [3 是 常用 的 斑点 噪声 滤波 算法 ， 其 方程 
Mclnerney 等 人 鲁 提 出 一 种 拓扑 自 适 应 主动 轮廓 模型 定义 为 
(T-Snake )， 随 着 轮廓 线 的 演化 ，T-Snake 模型 能 够 自 适 应 3) jegvI i (1) 
地 改变 其 拓扑 结构 ， 对 初始 位 置 的 要 求 较 低 。 备 令 奎 等 人 器 or 


提出 一 种 基于 正 交 T-Snake 模型 的 图 像 分 割 方法 ,Zhou 等 人 其 中 : 了 为 输入 图 像 ; 1 为 扩散 时 间 ; de 为 散 度 算 子 ，c(9) E 
09 提 出 基 结 合 全 局 与 局 部 分 布 信息 的 Snake 模型 的 医学 图 像 。” [0,1] 为 扩散 系数 方程 。 通常， 该 算法 对 与 中 心 像素 相 邻 的 四 
分 割 方法 ， 陈 能 成 等 人 0 提出 一 种 基于 T-Snake 模型 的 改进 ” 像素 进行 各 向 异性 扩散 计算 ， 其 缺点 在 于 图 像 特 征 边缘 常 被 
图 像 分 割 方法 。 虽 然 这 些 方 法 对 特定 的 图 像 分 割 效果 较 好 ， 过 度 平滑 ， 造 成 失真 。 基 于 窗口 的 改进 的 各 向 异性 扩散 算法 
是 还 是 存在 以 下 问题 ，a) 斑 点 噪声 对 于 分 割 方法 效果 的 影 (WSRAD) 5 用 子 窗口 中 的 像素 灰 度 均值 代替 四 邻 域 像素 
响 较 大 ，b) 边 界 连 续 性 差 、 边 界 特征 信息 弱 对 分 割 结果 影响 ， 灰 度 值 ， 其 缺点 在 于 方法 中 采用 的 中 值 滤波 无 法 保留 图 像 的 


卫 


民国 


较 大 ; c) 不 能 较 好 应 对 甲状 腺 超声 波 图 像 灰 度 不 均匀 、 边 缘 局 部 细节 信息 ， 亦 造成 失真 。 根 据 超 声波 图 像 的 斑点 噪声 的 

特征 信息 较 弱 、 轮 廓 连续 性 较 差 的 问题 。 分 布 特性 ， 本 文 在 WSRAD 滤波 方法 基础 上 对 子 窗口 进行 自 
针对 上 述 问 题 ， 本 文 提出 一 种 改进 的 T-Snake 模型 。 首 适应 加 权 中 值 滤波 ， 子 窗口 像素 (的 权 值 定义 为 

先 ， 结 合 自 适 应 加 权 中 值 滤波 算法 与 基于 窗口 的 各 向 异性 扩 W(i,)) =[@ -cdo? /ml] CO) 
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: 5 为 调节 因子 


等 ; 自 适应 渡 波 


口 像素 点 已 访 与 窗口 中 心 点 


的 距离 ， 5 为 子 
口中 心 像素 点 
取 中 心 点 权 值 吕 为 135，j 


(a) 初 始 区 


(a)Original image 


(c)WSRAD 滤波 
(c)WSRAD filter 


的 WSRAD 方法 很 好 的 保留 了 医 


口 灰 度 均值 ， 
为 取 整 运算 符 。 在 本 文 实验 中 
。 实 验 结果 如 图 


(b)SRAD 滤波 
(b)SRAD filter 


(gd) 本 文 方法 滤波 

(d)Proposed filter 
斑点 噪声 滤波 实验 结果 
Experimental results of speckle noise filtering 
Ph，SRAD 方法 产生 
降 块 效应 ，1 


效应 ，WSRAD 方法 避 
言 息 的 丢失 。 本 文 提出 


滤 除 了 斑点 噪声 ， 较 好 的 完成 了 超声 波 区 
改进 的 T-Snake 模型 
在 T-Snake 中 ， 其 


: 本 ,为 内 力 ， 其 表征 伍 
居 密 切 相关 ， 为 外 力 ， 表 全 


征 信息 ， 同 时 有 竟 
像 的 预 处 理 。 
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在 传统 的 Snake 模型 
= J a d+ BC GF ds 


: C(s) rie C Cs) 为 


像 梯 度 谋 ， 


es 性 等 图 
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线 OA AAA 
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模型 对 传 统 的 人 ,模型 


线 用 n 个 采 档 


点 表示 {x <isn,neNr*} 
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2 与 : 


4 为 da 与 di 的 加 权 平 均值 ; 


为 节点 的 相 邻 节 


进行 了 离散 化 处 理 。 
} ， 对 于 某 个 轮 


E, => 4 x -x -nb tbld; -di| 
i=1 
点 ; 4 为 相 邻 节点 之 间 的 


5 弹性 系数 ，2 为 刚 


E 质 密切 


线 的 一 阶 导数 ; 
4 弹性 系数 ; 
霄 能， 起 到 收缩 
刚性 势能 ， 起 到 PB 


B(s) 表 
线 的 作 
民 制 曲线 


性 系数 。 
常 月 
连续 
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有 T-Snake 模型 包括 脱 
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与 图 像 力 ， 对 于 图 像 中 清 
的 边界 轮 廊 检测 效果 较 好 ， 但 是 ， 在 甲状 腺 超声 波 
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a 加 | 
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小 方差 的 


为 膨胀 力 ， 使 得 轮廓 快速 收 龟 


其 中 ;Ejmoge 为 图 像 力 ， 


图 像 中 , 甲状 腺 轮廓 边缘 特征 信息 
本 文 在 传统 T-Snake 模型 基础 
区 域 能 量 项 [95]， 


, 连续 性 较 差 。 对 此 ， 


下 入 种 基于 局 部 区 域 最 
外 力 定义 为 


Es = YE + 
站 


Earion ) +nEre. gion ( 6) 


E 动 轮廓 向 


标 轮 廓 演化 逼近 Er ; 


项 ， 在 目标 轮廓 边界 特征 信 


又 域 能 量 项 保证 


权重 值 。 


区 


其 中 : < 为 


收敛 边界 的 ] 
像 力 Euee 定义 为 


ma 


标 轮廓 。 Ex 为 区 域 能 量 


息 较 弱 ， 连 续 性 较 差 的 情况 下 ， 
性 ; m 为 区 域 能 量 项 Exowio 的 


表示 [本 > 平 
膨 长 力 


其 中 : 


F(U(x,y)) 根据 轮廓 点 (x,y) 邻 切 


se =—C|VGG, y)* 17, y)| (7) 
V 表示 梯度 算 子 ; VG(x,y)*1(%,y) 
函数 与 输入 图 
定义 为 


Ennion = 9° F(T(x, y) n(x y) (8) 


图 像 统计 特性 决定 脱 


长 力 的 方向 nx,y) 为 轮廓 节点 (x,y) 处 的 法 向 量 大 小 ; 9 为 
自 适 应 的 膨胀 力 系 数 。 在 传统 方 当 9 通常 设 定 为 常数 ， 
为 了 达到 更 好 的 鲁 棒 性 与 收敛 速度 ， 本 文 定义 自 适 应 膨胀 力 
q=2—exp(-0—0o) (9a) 

0 a (9b) 

ac= 二 UND- (90) 


(xy)en 
0 为 轮廓 点 (xy) 的 3x3 邻 域 ， 轮 
区 域 ，5 为 邻 域内 最 大 灰 度 值 与 最 小 灰 度 值 
度 方差 ， 


.PH . 

内 部 区 域 与 外 部 
之 差 ; 5 为 邻 域内 灰 
线 位 于 灰 度 均匀 图 像 
有 特性 进行 演化 ; 轮 
2 之 间 ， 人 额外 的 脱 


区 域 能 量 项 Ercwo 


其 中 : 
万 内 部 区 域 灰 度 均值 ，& 为 轮廓 外 部 
部 区 域 ， 9 为 外 部 区 域 。 
从 数学 的 观点 看 ， 总 体能 
泛 函 ， 利 用 贪 禁 算法 可 求 得 总 
线 即 为 目标 轮廓 。 
3 ”实验 分 析 
3.1 实验 准备 
由 先 验 知识 可 知 ， 描 述 


廓 线 将 此 邻 域 分 割 为 


又 域 时 ， 

廓 位 于 灰 度 不 均匀 

长 力 驱 动 轮廓 加 速 离开 这 一 区 域 。 

定义 为 

ei ED- Up) | (10) 
LpeQ pe 


| 当 轮 廊 
此 时 轮廓 保持 原 
人 # 介 于 1 与 


ec 为 邻 域 9 内 点 p(x,y) 处 灰 度 值 ， 为 邻 域 在 轮 


区 域 灰 度 均值 ，9 为 内 


数 不 宜 


小 。 本 文 方法 中 ， 
系数 b=0.4; 


内 存 6 GB; 实验 数 ] 


模 


过 小 ， 刚 性 系数 不 宣读 
型 参数 i 


区 域 能 量 权重 系数 7=03 ; 
斯 平滑 函数 核 设 为 1.5 。 
实验 环境 为 : MATLAB 7.1, Lenovo X230i,CPU 2.4 GHz， 


昌 为 : 南通 市 第 一 人 民 医 院 提 供 的 20 套 


.最 小 时 的 轮 记 线 ， 此 上 


一 个 关于 轮廓 曲线 的 


人 廓 曲线 时 设 定 的 弹性 系 


或 能 量 权重 系数 不 宜 过 


单 性 系数 wa-=01， 刚 性 


图 像 力 系数 c=5 ; 高 
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甲状 腺 超声 波 图 像 ， 其 中 包括 了 原始 数据 集 与 专家 标注 集 ， 
对 数据 集 图 像 进 行 归 一 化 处 理 作为 实验 导入 数据 。 
3.2 评估 方法 
本 文采 用 超声 波 图 像 分 割 效 的 常用 评估 方法 :相对 习 
度 (ROD) 与 相对 差异 度 (RDD)。 两 者 公式 定义 为 

[CNCs| ICr NGC] 


uy 
路 


ROD = min pea x100% (11a) 
RDD = [CIC 00% (11b) a 
ICr (32) 甲状腺 肿瘤 超声 图 (b) 专 家 手动 分 割 
其 中 : CG 为 目标 轮廓 区 域 ,， 本 文中 以 专家 手动 分 割 结果 作为 (a)Thyroid ultrasonic image (b)Expert manual segmentation 
目标 区 域 ，C; 为 实验 收敛 轮廓 区 域 。ROD 表示 的 是 收敛 轮 


廊 曲 线 与 目标 轮廓 曲线 的 重 辣 程度 , RDD 表示 两 者 的 差异 程 
度 。 
3.3 ”实验 结果 
20 套 临床 数据 的 分 割 实验 结果 如 表 1、2 所 示 ， 

为 典型 的 分 割 实验 结果 图 示 。 
表 1 实验 结果 ROD 
Table 1 Experimental results-ROD 

最 大 ROD 最 小 ROD 平均 ROD 方差 


| 


2~4 


本 文 方法 94.3% 87.2% 91.3% 4.8% (c) 本 文 方法 (d) 正 交 T-Snake 
正 交 T-Snake 88.4% 76.4% 83.6% 9.2% (c)Proposed method (d)orthogonal T-Snake 
表 2 实验 结果 RDD 到 3 ”实验 数据 分 割 结果 II 
Table 2 Experimental results-RDD Fig.3 Results of segmentation II 
最 大 RDD 最 小 RDD 平均 RDD 方差 
本 文 方法 6.2% 3.6% 4.8% 1.6% 
正 交 T-Snake 12.7% 4.2% 8.6% 5.4% 


表 1 与 2 的 实验 结果 表明 ， 平 均 相 对 差异 度 小 于 5%， 
F 均 相对 重合 度 大 于 91%， 验 证 了 本 文 方法 的 可 行 性 。 
图 2~4 中 ， 本 文 方法 的 结果 与 专家 分 割 的 结果 吻合 度 较 
高 ， 较 好 地 处 理 了 甲状 腺 超声 波 图 像 灰 度 不 均匀 、 边 缘 特 征 
言 息 较 弱 的 问题 。 
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(a) 甲 状 腺 肿瘤 超声 图 (人 b) 专 家 手动 分 割 


(a)Thyroid ultrasonic image (b)Expert manual segmentation 


(a) 甲 状 腺 肿瘤 超 声 图 (b) 专 家 手动 分 割 


(a)Thyroid ultrasonic image (b)Expert manual segmentation 


(c) 本 文 方法 (d) 正 交 T-Snake 
(c)Proposed method (d)orthogonal T-Snake 
图 4 实验 数据 分 割 结果 王 

Fig.4 Results of segmentation II 
4 ”结束 语 
本 文 提 出 了 一 种 新 的 甲状 腺 超声 波 图 像 的 分 割 方法 。 在 


各 向 异性 滤波 方法 的 基础 上 ， 引 入 自 适 应 中 值 滤波 ， 有 效 地 
(c) 本 文 方法 (d) 正 交 T-Snake 滤 除 了 斑点 噪声 ;将 自 适 应 的 区 域 能 量 与 膨 外 力 引 入 
(c)Proposed method (d)orthogonal T-Snake T-Snake 模型 ， 实 现 了 对 弱 边 缘 特 征 信息 的 有 效 提 取 。 下 
图 2 ”实验 数据 分 割 结果 I 步 将 致力 于 将 本 研究 应 用 于 甲状 腺 超声 图 像 的 在 线 实 时 分 析 。 
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